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SECTION  1 

INTRODUCTION  AND  SUMMARY 


The  purpose  of  this  research  program  is  to  improve  our  understanding 
of  explosion  source  functions  and  source  coupling  through  theoretical  mcceiing 
of  underground  explosions,  and  through  laboratory  modeling  of  exclosicns  ana 
theoretical  simulations  of  the  laboratory  results.  Simulation  of  underground 
explosions  requires  a  detailed  knowledge  of  material  propert.es  and  const  tutive 
models  for  the  surrounding  medium  Some  of  these  properties  are  well 
constrained  by  observation,  while  others  are  much  less  certain.  In  this  reoort. 
we  summarize  our  current  understanding  of  the  constraints  on  mode!  ng  of 
underground  nuclear  explosions  in  granite. 

S-CUBED  has  a  history  of  very  successful  numerical  model. ng  of 
explosion  data  e  g  Rimer  and  Lie.  1982.  Rimer  ana  Cherry.  1S83:  Rimer  and 
Lie.  1984).  Near  field  data  from  the  explosion  PILEDRIVER.  for  example,  nas 
been  quite  accurately  reproduced  m  both  one-  and  two-dimensional  calculations 
(Rimer  et  al.,  1979;  Day  ef  a/.,  1983).  In  Section  2  of  this  report,  we  discuss  the 
models  that  were  used  in  these  simulations,  and  show  which  of  the  aspects  of 
the  models  are  required  by  the  data  and  which  parameters  are  less  well 
constrained.  We  then  scale  the  PILEDRIVER  'esults  to  correspond  to  the  y  eics 
of  HARDHAT  and  SHOAL  and  compare  the  results  with  near  field  wave'c'ms 
from  these  events 

In  Section  3  we  discuss  the  effect  of  d'latancy  ntrcduced  by  the  use  c‘ 
an  associated  how  rjie  in  simuiatons  of  explosions  in  g-amte  T^e  d.he-ence  m 
cavity  radu  between  the  Hoggar  expios-ons  and  NTS  exc  os'ons  could  ce 
explained  if  diiatancy  were  present  >n  the  Hoggar  exp. os  ons.  anc  absent  m  tne 
NTS  explosions.  However,  diiatancy  would  also  severely  change  the  gm-mc 
motion  The  ground  motion  in  the  diiatant  Hoggar  simulation  is  about  a  factor  of 
three  greater  than  the  ground  motion  .n  the  nondhatant  NTS  s.mihat  o° 
Diiatancy  also  has  the  effect  of  eliminating  expicsion  spectral  ^eaks  The  lame 
volume  increase  m  the  nonlinear  zone  near  the  expios.cn  makes  cwe’S-cut 
nearly  impossible 


S-CUBED  has  the  responsibility  *or  planning  anc  guidance  of  'abcra:ory 
experiments  now  in  progress  at  SRI  and  New  England  Research  Quasistatic 
tests  are  being  done  at  New  England  Research  to  determine  the  materia 
properties  of  the  Sierra  White  Granite  which  is  being  used  in  the  dynamic 
experiments  at  SRI.  In  Section  4.  we  report  on  the  current  status  of  tnese 
experiments. 
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SECTION  2 

CONSTRAINTS  ON  MODELING  OF  UNDERGROUND 
EXPLOSIONS  IN  GRANITE 


The  parameters  used  m  simulations  of  underground  ^uoear  exDiosions 
are  constrained  by  laboratory  tests  of  material  properties  and  by  direct 
observations  of  ground  motion  ‘rom  underground  explosions  In  this  secho~  we 
discuss  the  parameters  and  constitutive  models  that  were  used  in  S-CUBED 
granite  calculations  that  successfully  reproducec  near-field  wave'orms  and 
cavity  radii,  and  examine  the  effect  of  changing  these  parameters. 

Quasistatic  laboratory  tests  are  used  to  determine  the  fo!iowmg 
parameters:  density,  elastic  moduli,  pressure-volume  relation,  mater. al  strength, 
porosity,  and  water  content.  Of  these  quannt.es,  the  density,  moduli,  anc  P-V 
relation  seem  to  be  quite  consistent  with  in  situ  rock  properties  Porosity  and 
water  content  are  more  variable  m  situ  due  tc  the  presence  of  joints  and  are 
therefore  less  well  constrained  by  laboratory  data.  Laooratory  measurements  of 
strength  seem  to  be  very  inconsistent  with  n  situ  strength  as  interned  by 
modeling  of  underground  explosions. 

Direct  observations  of  underground  explosions  provide  seme  strong 
constraints  on  modeling.  In  particular,  cavity  radu  measurements,  near  f ■  e i d 
observations  of  velocity  and  diSDiacement.  and  m  cases  where  fu'l  ws.efc'ms 
are  net  available  (such  as  in  the  Hoggar  explosions',  peak  c  spiacemem  a-'d 
peak  velocity  measurements  an  can  be  compared  with  the  resuit$  of  nu^e'  cai 
simulations  of  undergrou.no  explosions  The  problem  wit"  usmg  p~ese 
observations  for  modeling  is  that  tne  data  quality  is  f'equent'y  inconsistent  a^d 
uncertain  Aiso.  unxnown  local  variations  in  material  properties  aoo  to  tne 
uncertainty  of  modeling  the  explosions  The  dynamic  experiments  tests 
described  m  Section  4  are  intended  to  be  simulations  of  uncerground  nuc  ear 
explosions  The  reason  for  performing  these  tests  is  to  try  to  'eoroduce  tne 
physics  of  underground  explosions  under  more  constrained  ccnoit  ons 


Complete  near  fieic  time  histones  are  available  from  three  Un  tec  States 
explosions  m  granite:  PILEDRIVER.  HARDHAT,  and  SHOAL  (Table  1)  We  nave 
digitized  the  near  field  data  from  these  events  for  comparison  with  our 
calculations  and  for  comparison  with  future  experimental  re  suits  Tne 
PILEDRIVER  data  was  collected  by  Sandia  Laboratory1 2  (Perret,  1968)  arc  by 
Stanford  Research  Institute^  (Hoffman  and  Sauer,  1969)  The  HARDHAT  data 
was  also  collected  by  Sandia  Corporation1  (Perret,  1963).  and  by  SRI  (Swift  and 
Eisler.  1962).  The  SHOAL  data  was  collected  by  Sandia  Corporation  (Wear:, 
1965).  Details  of  the  experiments  are  given  in  a  summary  by  Murphy  1 1973)  as 
well  as  in  these  reports. 

One-  and  two-dimensional  calculations  of  PILEDRIVER  performed  earlier 
at  S-CUBED  were  quite  successful  at  modeling  some  of  the  near  field  data  and 
reproducing  the  observed  cavity  radius  (Rimer  et  a!.,  1979  Day  et  at..  '963)  In 
addition  to  the  horizontal  data  shown  m  this  report,  there  is  a  substantial  amount 
of  vertical  data  available  from  PILEDRIVER.  and  this  data  was  reproduced  quite 
well  m  the  two-dimensional  calculations.  The  constitutive  models  usee  in  these 
calculations  are  discussed  m  detail  m  Rimer  et  al.  (1984).  In  Figures  i  anc  2.  we 
show  a  comparison  of  the  observed  near  field  velocity  and  displacement  from 
PILEDRIVER  compared  with  the  results  of  one-dimensional  calculation  No.  570. 
Toe  agreement  between  the  calculated  velocity  and  the  data  at  the  closest  two 
stations  is  exce'ient  At  the  more  distant  stations,  the  shape  of  the  velocity 
waveforms  is  reproduced  by  the  calculation  but  the  amplitude  of  the  waveforms 
is  overestimated  by  about  a  factor  of  twe  The  more  distant  SRI  stances  are  at  a 
o  f'erent  azimuth  tnan  me  closer  Sandia  stations,  so  tne  Observed  c  “erences 
could  oe  an  me  cation  of  azimuthal  var  aoons  c-f  tne  near  f  e:c  waveforms  As 
snown  m  Tabie  i .  tne  cavity  radii  from  tn  s  calculation  agree  very  wei  w  tn  tne 
observed  cavity  rad.i  (Heuze  "933)  for  ait  th-ee  events 

In  order  to  exam  re  me  amount  cf  variation  m  these  results,  we  scaled 
the  PILEDRIVER  calculation  to  the  vieid  of  the  explosions  HARDHAT  and  SHOAL 
and  plotted  the  resu'ts  together  with  near  f  eic  cbservat'pns  for  these  events 


1  Now  Sancia  Laboratories 

2  Now  SRI  International 


Velocity  (Meters/Second) 


Displacement  (Meters)  Displacement  (Meters) 


TABLE  1 


UNITED  STATES  EXPLOSIONS  IN  GRANITE 


Measured 

Cavity 

Explosion  Yield  (KT)  Radius  (m) 


Calculated  (570) 
Cavity 
Radius  (m) 


PILEDRIVER 

62.0 

40.1/44.5* 

42.5 

HARDHAT 

5.9 

19  4 

19.4 

SHOAL 

12.5 

26.8 

24.9 

*  44  5  meters  was  reported  in  Heuze  (1983).  40.1  meters 
was  obtained  by  a  drill  back  measurement  (Sterrett. 
1969) 
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HARDHAT  was  an  earlier  5  9  kiioton  explosion  located  very  c  ose  to  tne 
PILEDRIVER  site  m  the  same  medium  as  PILEDRIVER  The  HARDHAT 
comparison  is  shown  in  Figures  3.  4  and  5.  The  results  are  similar  to  the 
comparison  at  the  more  distant  PILEDRIVER  stations.  The  HARDHAT  peak 
velocity  amplitudes  are  overestimated  by  about  a  factor  of  two  at  the  Sanaia 
stations  and  by  about  a  factor  of  1.5  at  the  SRI  stations  however  the  calculated 
velocity  pulses  have  the  same  shape  as  the  data.  In  particular,  as  was  the  case 
with  PILEDRIVER.  the  observed  and  calculated  velocity  waveforms  are  broad 
and  have  large  negative  velocities  at  late  times.  The  HARDHAT  displacement 
data  at  the  SRI  stations  is  matched  very  well  by  the  calculation. 

SHOAL  was  a  12.5  kiioton  explosion  in  granite  located  away  from  tne 
Nevada  Test  Site,  so  the  material  properties  of  the  medium  are  not  identical  to 
those  of  PILEDRIVER  and  HARDHAT.  The  SHOAL  comparison  is  shown  in 
Figure  6.  The  results  are  similar  to  the  HARDHAT  comparison  Again  the  shape 
of  the  data  is  reproduced  fairly  well,  but  the  calculation  overestimates  the 
amplitude  of  the  velocity  pulses.  The  SHOAL  data  illustrates  one  of  the 
problems  with  such  comparisons  Note  that  all  three  observations  are  at 
approximately  the  same  range  (in  different  directions),  yet  the  observed  velocity 
waveforms  are  quite  different,  especially  at  station  PM-3  It  is  not  clear  whether 
the  data  is  actually  this  nomsotropic,  or  whether  there  is  some  problem  with  the 
instrumentation  that  recorded  the  waveforms  However,  as  we  remarked  earlier 
azimuthal  variations  were  observed  m  both  the  PILEDRIVER  data  and  the 
HARDHAT  data  as  well 

Peak  velocity  and  peak  displacement  data  from  explosions  in  granite  are 
also  available  from  tne  French  nuclear  explosions  m  the  Hocgar  (Heuze  1 983) 
The  Hoggar  data  and  tne  Deak  velocity  and  peak  displacement  cata  from  tne 
United  States  granite  tests  are  all  plotted  as  a  'unction  of  scaled  distance  tail 
data  has  been  scaled  tc  a  yield  of  one  kiioton)  m  Figures  7  and  8  The  Hcgcar 
near  field  data  is  quite  consistent  with  near  fie*d  data  from  the  Umted  States 
tests  Also  shewn  on  toese  figures  are  the  sca'ed  peak  velocity  and  peak 
displacement  curves  from  calculation  number  570  The  calculation  is  m  good 
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H8-RD  121m 


H8-RDI  121m 


HH-RD  1£4m 


H11-RDI  184m 


H12-RD  239m 


H12-RDI  239m 


4  Near  f ;  e 1  d  displacement  iron  HARDHAT  comparec  to 
d'Spiacement  from  calculation  nunoer  57C  scaled  to  5  S  kiiotons 
The  data  is  from  Perret  (iS63)  Tne  data  on  the  'eft  are  measured 
displacements,  and  the  oata  on  the  right  are  integrated  veioct'es 
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agreement  with  the  data  at  close  range,  however  as  we  saw  m  the  waveform 
comparisons  discussed  earlier,  the  peak  velocity  and  peak  displacement  cata 
are  lower  than  the  calculation  at  larger  ranges. 

In  the  calculation  (No.  570)  discussed  aoove.  aboratory  data  was  used 
to  constrain  the  elastic  moduli,  density,  and  the  initial  strength  of  the  rock 
(Cherry  and  Peterson,  1970).  The  failure  surface  used  m  the  caicuiat  cn  is 
shown  m  Figure  9.  The  strength  of  the  rock  at  P  equal  to  0  is  150  bars.  Toe 
strength  was  reduced  during  the  calculation  thro^gn  the  use  of  an  elective 
stress  law  (e.  g.  Rimer  ef  a/ ,  1984).  In  this  model,  the  medium  starts  out  a:  me 
laboratory  strength,  but  weakens  as  pores  are  crushed  by  the  explosion.  In 
porous  materials,  crushing  of  pores  can  have  a  very  strong  effect  on  source 
coupling,  however  the  porosity  of  granite  is  very  lew  (0.1  percent  maximum 
porosity  is  used  in  the  calculation),  and  the  main  effect  of  pore  crushing  is  to 
trigger  the  strength  reduction  of  the  effective  stress  law 

The  effect  of  this  model  is  to  severely  weaken  a  zone  around  the 
explosion  while  leaving  an  intact  high  strength  region  farther  out  As  shown  m 
Appendix  A.  a  weakened  core  inside  a  strong  rock  matrix  is  necessary  to 
generate  the  large  overshoots  and  negative  velocities  that  are  observed  m  the 
data.  It  is  well  known  that  if  laboratory  strengths  are  used  to  simulate  the  near 
field  data,  the  resulting  waveforms  are  much  too  narrow  and  lew  in  amphtjee 
(see  Figures  10  and  11).  and  the  resulting  cavity  radius  is  only  aoout  twc-tmrds 
of  the  observed  cavity  size.  In  Figures  12  and  13.  we  show  the  results  c*  a 
calculation  in  which  the  medium  :s  assumed  to  be  very  weak  ; 2 C 0  bar  ‘a;  ^re 
surface)  throughout  the  entire  calculation  The  result. ng  wavefems  are  a  fa  ' 
match  to  the  data  however  it  does  not  reproduce  the  necavve  veiocTes  at  me 
closer  stations,  and  the  calculated  cavity  radius  s  cverpred  ctec  oy  about  a 
factor  of  3- 2 

Either  the  elective  stress  law  or  some  otner  mechanism  for  se>e'e  y 
weakening  the  mater. ai  near  the  explosion  seems  to  be  requrec  to  match  cot" 
the  near  held  data  and  the  cavity  radii  It  is  not  clear  why  the  calculations  mate- 
some  of  tne  cata  very  well,  while  cverestimatmg  the  amplitudes  at  other  sta*  c-s 
Some  additional  numer.cal  mocei.ng  m  conjunction  w.th  the  exper -me-ta 
program  could  he'p  to  understanc  th  s  data 
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Near  field  velocity  data  from  PILEDRIVER  compared  with  cata 
from  a  calculation  using  a  strength  of  200  bars  throughout  the 
calculation 
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Figure  13.  Near  field  displacement  data  from  PlLEDRIVER  ccncarec  w  tn 
data  from  a  calculation  using  a  strength  of  200  Cars  throughout 
the  calculation. 
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SECTION  3 

THE  EFFECT  OF  DILATANCY  ON  GROUND  MOTION 
CALCULATIONS  IN  GRANITE 


Most  S-CUBED  ground  motion  calculations  have  used  a  racial  return 
nonassociated  flow  rule.  In  an  earlier  report  (Rimer  ef  a/..  1986).  we  showed  that 
using  an  associated  flow  rule  to  include  dilatancy  introduced  by  plastic  flow 
caused  significant  changes  in  both  cavity  radii  and  ground  motion  In  this 
section  we  address  the  questions  of  hew  an  associated  flew  rule  affects  the 
results  of  granite  calculations,  and  whether  dilatancy  could  explain  the 
Hoggar/PILEDRIVER  cavity  size  anomaly. 

Information  about  four  Hoggar  explosions  is  listed  in  Table  2.  When  the 
measured  cavity  radii  are  cube  root  scaled  to  a  yield  of  62  kdotons.  they  are  all 
about  30  meters,  much  less  than  the  measured  cavity  radii  for  the  three  United 
States  granite  shots  (Table  i)  which  all  exceed  40  meters  when  scaled  to  62  KT. 
Both  the  near  field  and  far  field  data  from  the  Hoggar  explosions,  however,  are 
comparable  in  amplitude  to  the  United  States  granite  data.  The  reasons  for  th'S 
difference  have  not  been  explained,  and  we  want  to  see  if  dilatancy  could  be 
responsible  for  this  anomaly. 

Dilatancy  cannot  be  added  directly  to  the  calculation  discussed  m  the 
last  section,  because  the  creation  of  a  large  volume  of  cracks  induced  by  p'ashc 
flow  is  clearly  incompatible  with  present  models  for  pore  crushing  ana  toe 
elective  stress  law.  To  examine  the  effect  of  dilatancy.  we  have  done  two  sets  of 
calculations  The  results  of  these  calculations  are  hstec  in  Tab'e  3  Figures 
showing  near  field  waveforms,  strain  paths,  and  reduced  d'Splacemert 
potentials  from  these  calculations  and  the  calculations  described  in  the  ast 
section  are  in  Appendix  B  First,  we  ran  the  calculation  with  the  high  laboratcry 
strength  (see  Figure  9)  with  an  associated  flow  rule.  The  effect  of  the  associated 
flow  -ule  is  largest  when  the  slope  of  the  failure  surface  is  large,  and  tnere  s 
difference  between  an  associated  flow  rule  and  a  radia1  return  ;iow  rule  when  the 
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TABLE  2 


HOGGAR  EXPLOSIONS 


Explosion 

Yield 

Measured 

Cavity 

Radius 

Cavity  Radius 
Scaled  to 

62  KT 

GEORGETTE 

13.0 

170 

28  6 

MICHELE 

3.6 

10  2 

263 

MONIQUE 

117.0 

390 

31.5 

CARMEN 

15  4 

200 

31  8 

22 

i 

:•> 

Si 

.vv  y-V  -y- 


1 1* * 
■*  *  a  m  ■ 


TABLE  3 


RESULTS  FROM  SIX  CALCULATIONS 


Calculation 

Cavity 

Radius 

(m) 

Elastic 

Radius 

(m) 

K 

(m3) 

^max 

(m3) 

No.  570 

42.5 

6  CO 

16000 

49000 

200  Bar  Strength 

62.6 

973 

53000 

57000 

Laboratory  Strength 
(Radial  Return) 

28.3 

385 

4700 

5000 

Laboratory  Strength 
(Associated  Flow  Rule) 

22.1 

580 

23000 

23000 

Slope  0  5  Failure  Surface 
(Radial  Return) 

41.3 

5*0 

15000 

15200 

Slope  0  5  Failure  Surface 
(Associated  Flow  Rule) 

30.5 

950 

58000 

58000 
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slope  of  the  failure  surface  >s  zero  Because  of  this,  dnatancy  is  small  '"■ear  tne 
explosion  because  stresses  are  very  high  and  failure  occurs  on  or  near  the  f:at 
part  of  the  failure  surface.  Dilatancy  increases  at  greater  distances  wnere 
stresses  are  lower,  and  will  also  occur  near  the  cavity  at  late  times  when  stresses 
are  reduced 

The  effect  of  dilatancy  is  to  slightly  increase  the  pulse  widths,  to  greatly 
increase  the  ground  motion  amplitude  especially  at  larger  ranges,  and  to  reduce 
the  cavity  size  In  Figures  14  and  15.  we  show  a  comparison  between  these 
results  and  the  near-field  PILEDRIVER  data.  The  waveforms  are  much  too 
narrow  and  are  a  poor  match  to  the  data.  The  cavity  radius  is  also  much  smaller 
than  the  observed  cavity  radius.  As  expected,  simply  adding  an  associated  flow 
rule  to  the  laboratory  data  does  not  produce  results  consistent  with  the  data 

In  the  second  set  of  calculations,  we  assumed  a  linear  failure  sjrface, 
and  varied  the  slope  and  intercept  of  the  curve,  searching  for  the  properties  that 
would  match  the  PILEDRIVER  cavity  radius  and  the  (scaled)  Hoggar  cavity 
radius  The  linear  failure  surface  leads  to  more  dilatancy  than  the  quadratic 
failure  surface  used  above  because  the  slope  of  the  linear  failure  surface  is 
steeper  at  higher  stresses  Each  calculation  was  run  twice,  once  using  a  radial 
return  flow  rule  and  once  using  an  associated  flow  rule  With  a  failure  surface 
slope  of  0  5.  an  intercept  at  zero  pressure  of  125  bars,  and  a  maximum  strength 
of  10  kilobars  the  calculation  with  the  radial  return  flow  rule  matches  the  cavity 
radius  of  PILEDRIVER.  and  the  calculation  with  the  associated  flow  rule  matches 
the  cavity  radius  of  the  Hoggar  explosions  The  ground  motion  in  the  diiatant 
ca:culation  (velocity  .  displacement,  and  RDP)  s  approximately  three  times  as 
arge  as  tne  grpjnd  motion  of  the  nondilatart  calculation 

Tne  fact  that  we  can  match  both  the  Hoggar  ana  PILEDRIVER  cav  ty 
radii  is  interesting,  but  it  is  net  clear  mat  it  has  any  'elation  to  reality.  Tne 
deserved  ground  motion  from  the  Hoggar  explosions  was  comparable  to  that 
from  PILEDRIVER  (Heuze.  1983).  while  the  calculations  show  that  cilatancy 
wejid  cause  much  larger  ground  motion  The  waveforms  m  ooth  calculations 
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Figure  14  Near  field  velocity  data  from  PILEDRIVER  comparec  with  cat 
from  a  calculation  using  the  laOoratory  strength  of  granite  arc  a 
associated  flow  rule 
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Figure  15 


are  a  poor  match  to  the  near  field  PILEDRIVER  ground  motion  (Figures  '6 
through  19).  It  is  unfortunate  that  there  are  no  near  held  ground  motion  records 
available  from  any  of  the  Hoggar  explosions  The  added  crack  volume  nduced 
by  the  associated  flew  rule  makes  overshoot  almost  impossible  As  a  result,  the 
large  negative  velocities  observed  near  PILEDRIVER  and  the  other  NTS  granite 
events  cannot  be  modeled  with  an  associated  flow  rule 

In  Table  3.  we  list  the  peak  RDP  value  and  the  fmai  RDP  value  for  each  of 
the  calculations  (figures  of  the  RDP  anc  RVP  for  eacn  calculation  are  in  Appencix 
B)  In  all  of  the  calculations  without  a  weaxeneo  core,  there  >s  \.ery  httie 
overshoot  (see  Appendix  A  for  a  theoretical  discussion  of  tne  relation  between 
the  weakened  core  and  overshoot).  The  effect  of  dilatancv  is  to  increase  tne 
amplitude  of  the  RDP  while  inhibiting  overshoot  even  further  Smce  the  increase 
in  amplitude  is  true  at  all  frequences,  an  explosion  that  generated  the  amount  of 
dilatancy  predicted  by  the  associated  flow  rule  would  cenerate  both  very  arge 
body  waves  and  very  large  surface  waves  compared  to  an  explosion  that  c>d  not 
produce  a  dilatant  region 
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Figure  16  Near  field  velocity  data  from  PILEDRIVER  compared  w>: h  ca:a 
from  a  calculation  using  3  fanure  surface  with  a  slope  of  0  5  ana  a 
radial  return  flow  rule 
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Figure  17  Near  field  displacement  data  from  PILEDRIVER  compared  v 
data  from  a  calculation  using  a  failure  surface  w>tn  a  ^icce  o* 
and  a  radial  return  flow  rule 
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Figure  18. 


Near  field  ve'OCity  oata  from  PILEDRIVER  comoared  witn  data 
from  a  calculation  using  a  failure  surface  witn  a  slope  of  0.5  ana 
an  associated  flow  rule 
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SECTION  4 

STATUS  OF  LABORATORY  EXPERIMENTS 


As  described  in  Section  2,  laboratory  quasistahc  material  properties 
measurements  on  granite  samples  show  a  much  larger  strength  than  that 
needed  to  simulate  the  PILEDRIVER  particle  velocity  measurements.  The 
strength  reduction  required  to  simulate  these  ground  motion  measurements  was 
assumed  to  be  due  to  pore  pressure  buildup  which  was  modeled  by  Rimer. 
et  ai,  1984  using  an  effective  stress  law  However,  other  mechanisms  for 
strength  reduction  are  possible  The  small  scale  explosive  experiments  being 
conducted  at  SRI  International  in  Sierra  White  Granite  (SWG)  are  intended  to 
isolate  the  physical  basis  for  the  near  field  ground  motion  measurements  from 
the  PILEDRIVER  event  in  NTS  granite.  As  described  by  Nagy  and  Florence, 
1986.  the  most  recent  senes  of  SRI  explosive  tests  did  show  a  factor  of  tnree 
increase  m  particle  velocity  pulse  widths  for  granite  test  samples  which  were 
extensively  gas  fractured  and  then  water  saturated  Less  fractured  samples  in 
earlier  tests  did  not  show  these  increases  in  pulse  widths  when  saturated,  nor 
did  gas  fractured  dry  samples.  The  more  fractured  laboratory  test  samples  are 
intended  to  more  closely  simulate  "in  a  scaled  sense"  the  in  situ  fractured  ana 
jointed  rock. 

As  described  in  Rimer,  etal..  1986.  major  problems  exist  m  interpreting 
the  experimental  results,  in  particular  the  small  cavity  radii  wmch  appear  to  oe 
inconsistent  with  the  velocity  measurements,  and  the  ack  of  rebound  impi.eo  oy 
some  of  the  panicle  velocity  measurements  It  is  also  dear  tKa:  the  two  halves  c* 
tne  cylindrical  test  samoles  had  separated  for  a  time  ai  owing  HE  detonaron 
or cducts  to  vent  along  tne  plane  of  the  part.de  velocity  gauges  ana  sweep  tne 
closest  m  gauges  out  of  their  emplacement  grooves  Recommencaticns  were 
made  which  shouid  improve  the  quality  of  the  expics  ve  tests  being  conducted 
this  year  at  SRI  Table  4  (from  Rimer.  et  al  .  1  586)  summarizes  the 
'eccmmended  explosive  test  senes  Numbered  tests  indicate  experiments 
completed  m  Fiscal  Year  1986 
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TABLE  4 

EXPLOSIVE  TEST  SERIES  IN  GAS  FRACTURED  SWG. 


Overburden  Pressure  Pore  Pressure 
Test  Type3  Pore  Condition  bars  (psi)  bars  f cs<; 


1  (473.  474)b 

Fully  Saturated 

136 

(2000) 

136 

■  2000 

2  (476)C 

Vacuum  Dried 

136 

(2000) 

0 

10) 

3d 

Fully  Saturated 

136 

(2000) 

3 

■50) 

4 

Fully  Saturated 

68 

(1000) 

3 

(50' 

5 

Vacuum  Dried 

68 

(1000) 

0 

(C) 

6 

Fully  Saturated 

136 

(2000) 

68 

i 1 CQC ! 

a.  Two  tests  are  recommended  for  each  test  type. 

b.  Fully  saturated  tests  were  actually  conducted  with  a 
vertical  confining  (clamping)  pressure  of  136  bars  but 
hydrostatic  overburden  and  pore  pressures  of  122  bars 
for  Test  473  and  109  bars  for  Test  474 

c  Test  476  did  not  have  complete  data  recovery. 

d  A  nominal  pore  pressure  of  a  few  bars  ;s  simcier 
experimentally  than  zero  pore  pressure 
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The  extensive  gas  fracturing  of  these  SWG  samples  clearly  aiter  tne.r 
material  properties.  Therefore  a  series  of  quasistatic  material  property  tests 
were  also  recommended  by  Rimer,  et  al.  (1986).  These  included  uniax.ai  stram 
load-unload,  triaxial  compression,  and  load-unload  tests  along  strain  patns 
relevant  to  the  SRI  experiments  The  material  properties  tests  are  being 
conducted  by  Dr.  Randolph  Martin  of  the  New  England  Research  Corporation 
On  August  18,  1966.  a  meeting  was  held  at  SRI.  with  S-CUBED  SRI  and  New 
England  Research  Corporation  represented,  to  coordinate  the  material 
properties  tests  with  the  dynamic  explosive  experiments  so  that  maximum  return 
is  achieved  from  the  limited  funds  available 

Granite  samples  fractured  at  SRI  have  already  been  received  by  New 
England  Research  and  Dr  Martin  anticipates  completion  of  the  first  phase  of  the 
material  properties  tests  by  the  end  of  October  SRI  has  recently  completed  the 
first  particle  velocity  experiment  of  this  new  test  series  This  experiment.  Test 
485.  was  a  repeat  of  dry  Test  476  with  deeper  grooves  for  all  the  particle  velocity 
gauge  wires,  including  comparisons  at  the  direct  ranges  between  these  grooves 
and  still  deeper  grooves  on  the  other  half  of  the  cylindrical  SWG  specimen  The 
sample  halves  again  separated  during  the  expenment  but  gauges  remained  in 
the  grooves  and  survived.  The  deepest  grooves  gave  the  cleanest  records 
Peak  displacements  and  cavity  size  agreed  well  with  Test  476  data  where 
available  and  filled  in  the  gaps  in  these  data 

Based  upon  this  successful  experiment,  SRI  will  next  execute  a  fjily 
saturated  experiment  with  very  low  pore  pressure.  Upon  completion  of  this 
experiment  and  the  fmst  phase  of  the  quasistatic  experiments,  another  meeting 
will  be  helc  at  SRI  with  S-CUBED.  SRI  anq  New  England  Research  Corporation 
representatives  present  to  analyze  resuits  before  proceeding  w  tn  the  remaining 
exDenments 
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APPENDIX  A 

SOME  ANALYTIC  RESULTS  RELEVANT  TO 
SHOCK  WEAKENING  IN  GRANITE 


The  explosion  source  function  (i.e..  the  moment  function  mft).  or 
equivalently,  the  reduced  displacement  potential  f(x))  can  be  related  to  source- 
region  mechanics  through  the  nonlinear  conservation  laws  and  the  constitutive 
equations.  The  latter  ought  ultimately  to  be  based  upon  laboratory  experiments 
inference  from  field  experiment,  micromechamcal  modeling,  or  all  of  these 

We  can,  of  course,  achieve  at  ieast  a  minimal  understanding  of  the 
seismic  source  without  actually  modeling  the  nonlinear  deformation  The 
cumulative  effect  of  source-region  deformation  is  to  exert  a  radial  stress  on  the 
spherical  surface  at  the  elastic  radius  r  As  a  first  approximation  one  can 
neglect  the  time  dependence  of  this  stress  pulse  and  treat  the  problem  in  which 
its  static  level  is  applied  instantaneously  This  is  the  classic  cavity  source 
problem,  whose  well-known  solution  is 


?(«) 


Pc(«] 


2 

u 


(A  1 ) 


where  P  (t)  is  the  pressure  applied  to  the  interior  cavty  wall,  starting  from  a 
quiescent  state  at  time  t  =  0  The  exterior  whoiesoace  is  isotroo  c. 
homogeneous,  linearly  elastic  The  exterior  wavespeeds  and  aens-t  es  S'e 
denoted  by  a,  p  and  p. 

This  source  spectrum  when  P  (t)  is  approximated  by  a  step  PqH(!)  is 
only  slightly  peaked  near  the  corner  frequency,  and  the  corresponding  far-fielc 
waveform  has  very  little  overshoot 

For  the  most  part,  dynamic  treatment  of  non, .near  ce'O'mation  within  the 
source  region  requires  numerical  sc'utions  to  the  govern, ng  camai  C'f,erentia. 
equations  However,  it  is  very  useful  to  have  a^a  ahe  a  wheal  re^u'*5  w-  c" 


aid  in  the  interpretation  of  numerical  solutions  In  particular,  it  would  be  useful  to 
have  some  insight  into  the  dynamic  behavior  engendered  by  a  weakened  shell 
of  rock  near  the  explosion  cavity,  abutted  by  stronger  rock  This  is  the  sort  of 
strength  contrast  we  would  expect  to  develop  dynamically  under  strain- 
weakening  or  shock-conditioning  conditions,  for  example 

For  the  purpose  of  developing  analytical  solutions,  these  circumstances 
can  only  be  represented  in  a  very  rough  manner.  By  analyzing  two  very  different 
idealizations,  however,  we  will  have  developed  some  interpretive  tools  with 
which  to  appreciate  more  rigorous  numerical  simulations. 

In  this  appendix,  we  first  investigate  weakening  effects  through  a  simple 
extension  of  the  linear  cavity  source  problem  Then,  nonlinear  effects,  both 
geometrical  and  material,  will  be  introduced  in  a  simplified  manner 

Neglecting  gravity  and  assuming  linear  elasticity  pertains,  we  consider  a 
spherical  bubble  of  radius  r1  in  the  center  of  a  fluid-filled  cavity  of  radius  r2  m  an 
elastic  wholespace  (Figure  A.1).  The  interior  and  exterior  wavespeeds  and 
densities  are  denoted  by  subscripts  1  and  2.  respectively 

The  bubble  is  pressurized  w.th  a  given  pressure  time  history  Pc(t).  The 
solution  for  the  displacement  s  then  has  the  frequency-domain  form 


s  =  -  U 
b  dr 


where 


^  -  i w r  a  -  iwr  a 

A  - - -  -  B  - - —  in  fluid. 


in  solid. 
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the  coefficients  A,  B  and  C  to  be  determined  from  tne  boundary  and  continuity 
conditions 


Trr  (r1>  =  -Pc 


T  rr  continuous  at 


s  continuous  at  r. 


The  radial  stress  in  terms  of  l  is 


frr  =  (X  -  2,)  |f  -  2X  f 


f, ,  ,  .  8^1  2X  Sj 

'  •  (x  ~  2l,)  ^2]  "  r  a? 


The  three  boundary/continuity  conditions  become 
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The  Cramer's  rule  solution  for  A  and  B  ts 


A  = 


-  b2  Pc 
alb2‘a2b1 


B  = 


g2  PC 

a1b2a2b1 


which  we  plug  into  Equation  (A.iO)  and  solve  for  C. 
simplifies  to 


C  = 
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The  expression  then 

(A. 20) 
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(A  23) 


42 


Equation  (A. 20)  gives  the  reduced  displacement  potential  wmch  we  nave 

/v 

heretofore  called  f  Thus,  s  is  given  by 


s  = 
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i  v 
a2rJ 


p  r  r  t 
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(r  * r2)  ' 1 
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2  q  (cos  6  -  a  sm  9) 


(A  24) 


In  the  high-frequency  limit,  Equation  (A. 24)  has  a  nice  mterpretat:on  We 
can  rewrite  it  as 
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The  high-frequency  limit  of  Equation  (A. 25).  for  the  case  of  P  equal  to  a  step 

V 

function  of  amplitude  P0  is 
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where 


r  =  p2q2  ;  p^a^ 
p2a2  ~  p^a^ 


(A .27) 
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Equations  (A. 27).  (A  28)  and  (A. 29).  respectively,  are  the  normal  incidence  plane 
wave  reflection  coefficient  (R)  and  transmission  coefficient  (T)  at  the  fluid-soiid 
boundary,  and  the  one-way  phase  delay  [6)  across  the  fluid  shell  Thus, 
Equation  (A. 26)  has  an  obvious  interpretation.  At  sufficiently  high  frequency,  the 
radius  of  curvature  of  all  boundaries  eventually  becomes  negligible  comoared 
with  the  wavelength  The  series  terms  then  represent  multiple  reflections  at  the 
fluid-solid  interface  The  first  motion  of  each  reflection  is  given  by  the  inverse 
transform  of  the  corresponding  series  term  (Abel's  theorem)  giving  the  following 
first  motion  senes,  most  easily  expressed  m  terms  of  velocity: 


u(r-'>  ■  yvi  T  £  (•>"  R"  h  ft  ■ 

Ha1r  n=0 


1  [  r2 "  rJ 


(A  30! 


where  H  is  the  unit  step  function. 

Thus,  the  fluid  shell  gives  rise  to  repeated  overshoot  and  rebound  m  the 
motion,  as  well  as  the  potential,  at  r^.  This  tends  to  put  a  stronger  peak  into  the 
far-field  displacement  spectrjm.  In  a  real  material,  there  would  be  some 
residual  strength  after  strain-weakening  (or  shock  conditioning,  or  wnatever 
weakening  mechan  sm  is  postulated  to  operate),  and  this  would  rapidly  damp 
the  oscillation  of  the  source  region,  probably  after  a  single  cycle  of  outward 
motion  and  rebound  Nonlinear  numerical  simulations  with  shock  weakening 
(effective  stress  law.  for  example)  indeed  exhibit  this  enhanced  rebound  and 
spectral  peaking 

Figure  A. 2  shows  some  complete  time  histories  for  tne  reduced 
d'splacement  potential  for  the  fluid-shell  problem,  in  which  the  cversncot  anc 
rebound  are  evident  The  extent  of  overshoot  seems  to  increase  with  the  extent 
of  the  fluid  shell  region 

The  preceding  solution  to  the  linearized  momentum  equation  with  elastic 
stress-strain  relation  is  useful  in  its  own  right,  m  that  it  demonstrates  that  strain 
weakening  near  the  source,  represented  here  is  an  approximate  way  oy 
mtroducing  a  fluid,  (zero-strength)  interior  snell.  can  lead  to  reverberative 
response.  As  noted  above,  it  seems  likely  that  if  the  interior  region  were  a  soft, 
plastic  shell  instead  of  a  fluid,  one  might  still  see  the  outgoing  pulse  (first  term  in 


Equation  A  30)  with  somewhat  diminished  amplitude,  the  f  rst  'ebound  [second 
term  m  Equation  A  30)  also  diminished  in  amplitude,  anc  then  locking  m  c*  tne 
cavity  displacement  due  to  the  strength  of  the  plastic  shell  This  would  be 
analogous  to  the  oscillation  followed  by  locKing  of  a  narmomc  oscillator  to  which 
a  Coulomb  slider  has  been  added.  If  this  qualitative  description  is  correct,  then 
strain  weakening  should  lead  to  a  source  function  with  significant  overshoot  To 
ma*e  this  notion  more  quantitative,  we  seek  the  help  of  an  analytic  solut  on 
Tms  calculation  is  a  use^jl  supplement  to  the  last  section  as  it  incorporates  both 
geometrical  and  material  nonlinearity  The  price  paid  is  that  the  source  region  is 
treated  as  incompressible,  which  leads  to  unphysical  behavior  of  the  first  motion 
Furthermore,  we  manage  only  to  reduce  the  problem  to  that  of  solving  a  system 
of  three  coupled  first-order  ordinary  differential  equations.  These  are  easily 
solved  numerically,  however. 

Figure  A  3  shows  a  spherical  shell  of  interior  radius  r  (t)  and  exterior 
radius  r2(t).  both  of  which  are  time  varying.  The  left  hand  side  is  the  initial 
configuration,  at  t  =  0.  the  right  hand  side  is  the  configuration  at  some  later  time 
t.  The  shell  consists  of  a  rigid-perfectly  plastic,  incompressible  solid  m  welded 
contact  with  a  linearly  elastic  whoiespace  at  r2  The  shell  has  density  f  1  and 
yield  strength  Y.  the  whoiespace  has  density.  P  wavespeed  and  S  wavespeed 
equal  to  p2.  a 2  and  /?2  respectively 

We  wisn  to  fmo  the  seismic  radiation  when  the  shell  is  leaded  at  time 
cero  by  a  pressure  Pc  which  may  .n  turn  depend  noninear'y  _pcn  the  motion. 

The  radial  coordinate  m  the  initial  unceformed  conf  curaton  will  be  our 
re;erence  cr  Lagrangian.  coordinate,  denoted  H  The  correspond. ng  spatial 
coordinate  will  be  denoted  r  =  r(R.t)  sc  that r  R  O'  =  R 

From  the  .ncqmpress.btlity  condition  we  have 

p(R.t)  -  p  R.0)  for  R1  <  R  <  R2  (A  31) 

anc  the  continuity  equation  under  spherical  symmetry  combmed  with  Equation 


The  plastic  yield  condition  is 


'W!Y 


'A  38: 


Because  of  the  rigid-plastic,  incompressible  conditions  there  can  be  no 
motion  unless  the  whole  plastic  shell  is  at  the  yield  surface,  then 


P  ^-§  =  -  7  sgn  (u)  Y  . 

dr 


(A  39) 


the  solution  to  which  we  seek,  with  the  boundary  conditions 


Vv 
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Trr(r2) 
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The  outer  boundary  stress  T 2  is  to  be  determined  by  requiring  continuity 
between  the  interior  plastic  solution  and  the  exterior  elastic  solution,  which  is  just 
the  cavity  source  problem  with  Pc  =  -T ^  Substituting  Pc  =  - T 2  into  Equation 
(A.1),  and  making  the  operational  interpretation  i u  =  d  dt.  leads  to  the  following 
ordinary  differential  equation  for  s2  =  s(r2.t): 


(A  42; 


If  we  can  express  the  solution  to  Equation  (A. 39)  in  terms  of  the  unknown 
boundary  value  T2.  then  that  expression,  combined  with  Equation  (A. 42).  wT 
constitute  an  ordinary  differential  equation  system  which  we  can  solve  for  s2  anc 
Tj.  Then  the  source  function  f  in  terms  of  s2  is  given  by  the  ccovoiuVcn. 


>  V  V  V 


To  solve  Equation  (A  39).  we  integrate  with  respect  to  r 


Tr  r  ( r )  =  -  Pc  -  2Y  sgn  (u)  In  f ~  -  p 


u(  r ' )  dr 


(A  44; 


Next,  evaluate  the  final  term  m  Equation  (A. 44)  in  terms  of  s2-  The  velocity  is 


u(r,t)  =|f 


R  =  R(r) 


(A  45) 


where  R(r)  is  the  inverse  function  of  r(R).  From  Equation  (A  34).  we  get 


=  u(r . t)  =  s. 
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(A  46) 


The  acceleration  is 


d  s  du  3u 
a,2  -dt  -  3 1 
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(A  47 


Plugging  Equation  (A  47)  into  Equation  (A. 44)  and  carrying  out  the  integration 
leads  to  an  expression  for  Trr(r,t),  which  we  evaluate  at  r  =  r2. 
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T2(t)  =  -  Pc-2Vsgn  (s2)  ^  «r 


2  2"1 
s*  - 1  s*j 


(A  48i 


where  5  r2  -  rr 


We  could  substitute  Equation  (A. 48)  into  Equation  (A. 42)  to  obtain  a  third-order 
ordinary  differential  equation  for  s2(t).  This  equation  is  comp  eted  by  specifying 
P  (r^)  via  an  equation  of  state,  e  g.,  the  constant-gamma  gas  approximation 


fR  ^7 

p  =  3XUjW 
C  4 r  R^3  lri. 


(A. 49) 


where  W  is  the  explosion  energy.  Equation  (A  48)  is  nonlinear  in  s2  and 
therefore  Equation  (A. 42)  generally  will  require  a  numerical  method  for  <ts 
solution. 

A  more  convenient  form  for  numerical  treatment  ensues  i*  we  define  the 
new  dependent  variable  r, 


t  =  s2  - 


4  P 2  a2 

Then  we  define  the  vector  components  V  by 


Vi2t 
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With  these  definitions.  Equations  (A. 42)  and  (A  48)  are  equivalent  to  the  following 
system  of  first  order  ordinary  differential  equations 


dV 

df-E  (¥•*) 
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(A. 53c) 


Recall  that  ^  is  a  known  function  of  s2,  given  by  Equation  (A  34)  toy  virtue  of  the 
incompressibility  condition),  and  that  P„  r  )  is  given  by  Equation  (A. 49). 
Equation  (A  52)  is  easily  solved  using  a  fourtn-orcer  Runge-Kutta  method,  with 
due  caution  regarding  tne  sign  of  V2  in  the  fourth  term  of  Equation  A  53b) 

Because  of  the  mcompressio  iity  condition  Equation  (A. 52)  is  unphysical 
m  its  representation  of  the  first  motion,  and  probably  generally  smooths  out 
high-frequency  features  which  would  appear  if  compressibility  were  permitted. 
This  smooth.ng  is  evident  in  numerical  results  shown  in  Figure  A. 4.  In  all  cases 
shown,  the  explosion  energy  is  50  KT,  emplaced  in  a  cavity  of  initial  radius  10  m 
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density  and  exterior  wavespeeds  are.  respectively.  2500  kg,m  ,  5000  msec,  and 
3000  m  sec.  Results  are  given  in  terms  of  the  reduced  displacement  potential, 
via  Equation  (A. 43)  The  first  motion  is  quadratic  in  time  rather  than  linear,  an 
artifact  of  the  incompressibility  condition.  When  the  shell  is  weak,  the 
anticipated  overshoot  appears,  and  damps  out  rapidly,  as  expected.  When  the 
shell  is  weak  enough  that  several  half-cycles  appear,  the  oscillation  amplitude 
diminishes  by  an  approximately  constant  Decrement  at  each  half  cycle  This 
behavior  roughly  mimics  that  of  a  harmonic  oscillator  m  parallel  with  a  Coulomb 
slider.  Thus,  m  spite  of  some  severe  oversimplification,  this  solution  verifies  that 
shock  conditioning  is  responsible  for  introducing  or  enhancing  overshoot  m  the 
source  function 


APPENDIX  B 

STRAIN  PATHS,  TIME  HISTORIES  AND  REDUCED  DISPLACEMENT 
POTENTIALS  FOR  GRANITE  CALCULATIONS 


"VP‘'""3>  VKodty  (ft^nc) 
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Figure  B  i.  Reduced  displacement  potential,  reduced  velocity  potential,  and 
displacement  and  velocity  at  1 2 2 5  meters  for  PILEDRIVER 
Calcu.ation  Number  57C 
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Figure  B  2 


Displacement,  velocity  stress  and  strain  time  h, stones  at  204 
meters  for  PILEDRIVER  Calculation  Numoer  570 
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Figure  B  4  Reduced  displacement  potential,  reduced  velocity  potential,  and 
displacement  and  velocity  at  1225  meters  for  weak  ca  cu  at  on 
strength  =  200  Pars) 
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Figure  B  5  Displacement,  velocity,  stress  and  strain  time  histories  at  204 
meters  for  weak  calculation  (strength  =  200  Pars) 


60 


volumetric  strains  (%) 


mean  stress  (KB) 


Figure  B  6 


Strain  paths  at  204  meters  for  weak  calculation  ^strength 
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Figure  B  7 


Reduced  displacement  potential  reduced  velocity  potential.  ana 
displacement  and  velocity  at  1 225  meters  for  nigh  iaooratory 
strength  calculation  -  radial  return  flow  r_:e 
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Figure  B  8  Disolacement.  velocity,  stress  and  strain  t.me  histories  at  234 
meters  for  high  laboratory  strength  calculation  -  rad. a l  return 
flow  rule 
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Reduced  displacement  potential,  reduced  velocity  potert  ai 
and  displacement  and  velocity  at  i  225  meters  for  mgr 
laboratory  strength  calculation  -  associated  flow  rje 
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Strain  paths  at  204  meters  for  high  laboratory  strength 
calculation  -  associated  flow  rule 
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Reduced  displacement  potential  reduced  velocity  potential, 
and  displacement  and  velocity  at  1225  meters  for  calculation 
with  a  linear  failure  surface  win  a  slope  of  0  5  -  radial  re’urn 
flow  rule 
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Figure  B  14  Displacement,  velocity,  stress  and  strain  time  histories  at  2C4 
meters  for  calculation  with  a  linear  failure  sur'ace  with  a  sicoe 
of  0  5  -  radial  return  flow  rule 
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Figure  B  16 


Reduced  displacement  potential,  reduced  velocity  potential 
and  displacement  and  velocity  at  1225  meters  for  ca'cu'ation 
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